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ABSTRACT 

We report results from four recent observations of the ultracompact low-mass X-ray binary pulsar 
4U 1626—67. All the observations obtained high-resolution X-ray spectra of the system, two from the 
Chandra X-ray Observatory using the High Energy Transmission Grating Spectrometer, and two from 
the XMM-Newton Observatory using the Reflection Grating Spectrometer as well as the EPIC PN and 
MOS. These data allow us to study in detail the prominent Ne and O emission line complexes which 
make 4U 1626—67 unique among LMXBs. The observations were spaced over a period of 3 years for a 
total observing time of 238 ks, allowing us to monitor the line regions as well as the overall source flux, 
continuum spectrum, and timing properties. The structure of the emission lines and the ratios of the 
components of the helium-like Ne IX and O VII triplets support the hypothesis that they are formed 
in the high-density environment of the accretion disk. We do not find any significant changes in the 
line widths or ratios over this time period, though we note that the line equivalent widths decrease. 
Using the most recent calibration products, we are able to place constraints on the strengths of the 
Ne K, Fe L, and O K photoelectric absorption edges. In contrast to our earlier analysis, the data do 
not require an overabundance of Ne or O in the system relative to the expected ISM values. We find 
that the pulsar is still spinning down, though the rate of decrease is greater than what was predicted 
by the ephemeris derived from previous Comp ton/BATSE monitoring. We also note that the pulse 
profile has changed significantly from what was found prior to the torque reversal in 1990, suggesting 
that this event may be linked to a change in the geometry of the accretion column. The flux of 4U 
1626—67 continues to decrease, in keeping with the trend of the last 30 years over which it has been 
observed. Taking into consideration current theory on disk stability, we expect that 4U 1626—67 will 
enter a period of quiescence in 2-15 years. 

Subject headings: binaries: close — stars: neutron — pulsars: individual (4U 1626-67) — X-rays: binaries 



1. INTRODUCTION 

The 7.7 s X-ray pulsar 4U 1626—67 was first discovere d 
by Vhuru (jGiacconi et al.lll972l: iRappaport et al.lll977D . 
and remains the only k nown high- field (i? w 4 x 10^^ G, 
IPravdo et al. 1979; Co burn et 311120021 ) pulsar in an ul- 
tracompact low-mass X-ray binary (LMXB). This unique 
pairing of an apparently young neutron star with a very 
low-mass companion could indicate that the neutron star 
was originally a white dwarf whose mass, as a result of 
accre tion, exceeded the Chandrasekhar limit (i Joss et ahl 
119781 though see Verbunt et al. 1990). Ahhough orbital 
motion has never been detected in X-ray data, pulsed 
optical emissioi i reprocessed on the sur face of the sec- 
ondary allowed [Middleditch et al.l (|1981h to infer an or- 
bital period of 42 min, which was later confirmed by 
IChakrabartyl ()1998l ) . It is therefore a member of the class 
of objects known as "ultracompact" binaries (Porb < 80 

^ Also Department of Physics, Massachusetts Institute of Tech- 
nology, Cambridge, MA 02139 



min) which must have hydrogen-depleted secondaries to 
reach such s hort periods (Paczynski & Sicnkiewicz 1981]; 
iNelson et al.lll986t ) . It has an extremely small mass func- 
tion of / < 1.3 X 10~^ M0, which corr esponds to a 
secon dary mass of 0.04 Mq for i — 18° (|Levine et al.l 
mil). A very low-mass secondary would account for the 
faint {V K, 17.5) op tical counterpart and the high optical 
pulsed fraction (McCh ntock et al.lll977l . IT980f) . 

Initially, 4U 1626—67 was observed to be spinning up 
with a characteristic timescale P/P k, 5000 yr, but 
in 1990, this trend reversed and the neutron star be- 
gan to spin down on approximately the s ame t imescale 
(IWilson et all IT993t iChakrabartv eTall [1993). The 
torque reversal was abrupt, although the decrease in 
bolometric X-ray flux has been gradual and continuous 
over the past k, 30 yr. Assuming that the X-ray flux 
is a good proxy for the accretion rate onto the neu- 
tron star, these observational facts cannot be reconciled 
with current accretion disk-neutron star interaction the- 
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ory. However, the fact that the pulsar underwent torque 
reversal implies that the radius at which the accretion 
disk is truncated by the neutron star's magnetic field 
is close to the corotation radius. With this assumption 
and a neutron star mass of 1.4 M©, the measured P 
during spin-up can constrain the accretion rate to be 
^, 2x 10~^° Mf7) yr~\ whi ch implies a distance of > 3 kpc 
(jChakrabarty et aTlllOQTl ). Alternatively, measurements 
of the optical and X-ray fluxes (assuming the albedo of 
the disk is > 0.9) give a distance ran ge to 4U 1626—67 of 
5 < < 13 kpc (|Chakrabartvl[l998h . 

Perhaps the most unique characteristic o f 4U 
1626—67 is its X-ray spectrum. lAngelini et alj ()1995 f) 
first reported the presence of Ne and O emission lines in 
an ASCA spectrum, and they found similar features in 
an Einstein observation performed prior to the torque re- 
versal. Using a subsequent observation by Chandra, we 
discovered double-peaked line st ructure indicat ive of for- 
mation in an accretion disk (Sch ulz et al.ll200 "ll. hereafter 
SOI). Emission line features such as these are not seen in 
any other LMXB systems, and suggest that the donor is 
particularly rich in elements resulting from later stages of 
nuclear burning — perhaps a C-O-Ne or 0-Ne-Mg white 
dwarf. A UV spectrum obtained with ffST/STIS re- 
vealed bot h emission and ab sorption features from C, 
O and Si (jHomer et al.l[200^ . The C absorption lines 
are stronger than expected from a purely ISM contribu- 
tion, suggesting some local contribution, and the O V 
line has a pronounced double-peaked profile, indicating 
an accretion-disk origin similar to the X-ray lines. The 
spectrum is missing common N and He lines typically 
seen in UV spectra from high-excitation systems, likely 
because 4U 1626—67 is lacking in th ese elements. A high - 
S/N optical spectrum obtained bv IWerner et al.l (|2006D 
using the VLT confirms the lack of He in the system, and 
finds H lacking as well. The optical spectrum is domi- 
nated by C and O emission lines, but does not appear to 
contain any Ne lines. Given the current understanding 
of the model Ne atom, they find that Ne is present in 4U 
1626—67 at < 10% by mass; or, equivalently, that Ne/0 
^ 0.2. 

We present a series of four high-resolution X-ray spec- 
tra spanning three years, two obtained with the High En- 
ergy Transmission Gratings (HETGS) aboard the Chan- 
dra X-ray Observatory, and two from the Reflection 
Grating Spectrometers (RGS) aboard XMM-Newton. 
We previously presented analysis of the first Chandra ob- 
servation in SOI, but this re-analysis uses improved cal- 
ibration products and software. We describe the obser- 
vations and data reduction procedures in Section [21 and 
our spectral analysis in Section [31 We present timing 
analysis is Section [4l Finally, we discuss the implications 
of our observations in Section [5l 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Chandra 

4U 1626-67 has been observed twice with the 
HETGS on board th e Chandra X-Ray Observatory 
(jCanizares et al.|[2005D . first on 2000 September 16 and 
again on 2003 June 3. The HETGS comprises two 
sets of transmission gratings: the Medium Energy Grat- 
ings (MEGs), with a spectral resolution of AA — 0.023 
A FWHM and a range of 2.5-31 A (0.4-5.0 keV), and the 



High Energy Gratings (HEGs), which have AA = 0.012 
A FWHM and a range of 1.2-15 A (0.8-10 keV). See 
Table [H for a summary of the observations. 

We obtained "level 1" event lists from the Chandra 
data archive^, and reprocessed the data using the latest 
available version of the CIAO software (v. 3. 3) and the 
calibration database (CALDB v. 3. 2.1). This processing 
applied the gain, time-dependent gain, and charge trans- 
fer inefhciency (CTI) corrections. The CTI corrections 
improve the computation of event grade (used to filter 
out likely non-source events), and the gain products im- 
prove the channel to energy mapping of events, which al- 
lows for more precise order-sorting of the dispersed spec- 
tra. We created grating responses that include the time- 
dependent effects of a contaminant present on the ACIS 
optical blocking filter ([Marshall et al.l[2004[ ). Correcting 
for this contamination is necessary for precise spectral 
analysis, since it affects both spectral shape and normal- 
ization. 

Neither Chandra observation contained any apprecia- 
ble background flares. After reprocessing and filtering 
out bad pixels, we applied the standard grade filtering 
(retaining grades 0, 2, 3, 4 and 6), and used the tgdetect 
tool to determine the zeroth-order source position. The 
high source count rate caused the zeroth-order images to 
be quite piled up, but this did not hinder our ability to 
determine the source position. This position was used 
to sort the dispersed spectra into the proper orders, re- 
sulting in 62555 (101799) background-subtracted events 
in the first-order MEG spectrum and 35343 (59706) in 
the first-order HEG spectrum for ObsID 104 (3504). The 
dispersed spectra are not affected by pileup. 

2.2. XMM-Newton 

4U 1626-67 has been observed four times with XMM- 
Newton, but only two of these observations contain a 
significant amount of science data (ObsIDs 0111070201, 
performed 2001 August 24, and 0152620101, performed 
2003 August 20; see Table [TJ. Here, we analyze data 
from the three EPIC cameras (two MOS detectors and 
the PN), as well as the two reflection grating spectrom- 
eters (RGSl and RGS2). The MOS (PN) detectors have 
nominal bandpasses of 0.15-12 (0.15-15) keV, and spec- 
tral resolutions of 70 (80) eV at 1 keV. The RGSs have 
bandpasses of 0.35-2.5 keV and first-order spectral res- 
olutions of 0.04 A FWHM. Both observations were per- 
formed after the failure of RGSl's CCD7 and RGS2's 
CCD4, resulting in gaps in the dispersed spectra which 
were filtered out for analysis. 

The XMM-Newton science products were obtained 
from the HEASARC^ archive. Both observations were 
affected by strong background flares in the MOS and PN 
detectors, which we filtered out prior to analysis, result- 
ing in the loss of w 5 ks of MOS data and « 2 ks of PN 
data for ObsID 0111070201, and « 25 ks of MOS and PN 
data and « 12 ks of RGS data for ObsID 0152620101. 

During ObsID 0111070201, the MOSl detector was in 
timing mode, so was not easily susceptible to pileup. 
However, the source core is somewhat piled in the M0S2 
detector, and was excised prior to analysis. During Ob- 
sID 0152620101, the M0S2 detector was in timing mode, 

^ http:/ /cda. harvard.edu/chaser/mainEntry.do 

^ http:/ /heasarc.gsfc.nasa.gov/db-perI/W3Browse/w3browse.pl 
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TABLE 1 
X-Ray Observations 



Observatory Observation start Observation ID Duration 



Chandra 2000 September 16 14:57 UT 104 40 ks 

XMM-Newton 2001 August 24 02:57 UT 0111070201 17 ks 

Chandra 2003 June 3 02:30 UT 3504 97 ks 

XMM-Newton 2003 August 20 05:55 UT 0152620101 84 ks 



and the MOSl detector shows mild evidence of pileup, so 
we again removed the central portion of source data prior 
to analysis. The standard filters were applied, FLAG 
== and PATTERN <^ 12 (4), to the MOS (PN) 
data. The RGS data were reprocessed (using XMM- 
SAS v.6.5.0) with refined source coordinates, resulting 
in updated source and background spectra as well as re- 
sponses. Background spectra were created for the EPIC 
data using off-source regions, and scaled appropriately to 
match the source spectra. After processing, there were 
a total of 1.1 X 10^ (1.7 x 10^) background-subtracted 
counts in the MOSl data, 7.3 x IQ-* (4.5 x 10^) in the 
M0S2 data, 3.1 x 10^ (1.04 x 10*^) in the PN data, 8366 
(17289) in the RGSl spectrum, and 9951 (21076) in the 
RGS2 spectrum for ObsID 0111070201 (0152620101). 

3. SPECTRAL ANALYSIS 

We used the ISIS software package'* (vl.3.3) for spec- 
tral fitting. All source and background data, as well as 
response files, were read into ISIS, and normalized back- 
ground counts were subtracted prior to fitting. We com- 
bined the ±1 spectral orders of the HETGS data, and 
used the —1 order of the RGS data. We employed Cash 
statistics throughout our analysis, and fit the data from 
different instruments simultaneously. All the quoted er- 
rors are 90% confidence limits unless otherwise noted. 



3.1. Continuum Fitting 

For continuum fitting, we excluded the known emis- 
sion line regions, and used the wavelength/energy ranges 
1.2-17.0 (1.8-26.0) A for the HEG (MEG) grating spec- 
tra, 7.0-34.0 A for the RGS grating spectra, and 1.1-9.9 
(1.1-8.3) keV for the MOS (PN) spectra. The lower en- 
ergy ranges of the MOS and PN spectra were excluded 
because they contain unresolved Ne and O emission lines 
which artificially increase the overall continuum level. 
We also excluded regions in the RGS data that contain 
bad columns or fall across non-functioning CCD chips. 
During the second XMM-Newton observation (ObsID 
0152620101), M0S2 was in timing mode, and its spectral 
calibration did not match that of the other instruments, 
so these data were omitted from the continuum fit. 

We grouped the MOS data to have a minimum of 50 
counts and the PN data to have a minimum of 500 counts 
per bin, and the RGS data to contain 6 channels per bin 
(resulting in ~ 0.04 A wavelength bins, comparable to 
the instrumental resolution). We grouped the HEG and 
MEG spectra to contain 4 channels and a minimum of 10 
counts per bin (resulting in > 0.01 [0.02] A wavelength 

4 http://space.mit.edu/CXC/ISIS/ 
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Fig. 1. — The combined Chandra (top panel) and XMM- 
Newton (bottom panel) spectra with representative continuum 
models comprised of absorbed blackbody plus pov^er-law emission. 
The blackbody (dot-dashed line) and power-law (dashed line) con- 
tributions are also plotted. 



bins for the HEG [MEG] data, again comparable to the 
instrumental resolution). 

Although a simple power-law model gives a reasonable 
fit {x1 — 1.2), the residuals suggest the presence of an 
additional spectral component. This has previously been 
modeled as a second power-law or a blackbody compo- 
nent (SOI). An F-test indicates that the addition of a 
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TABLE 2 
Continuum Spectral Fits 







Power-law 


Blackbody 






Observation 


{102lcm-2) 


Norm-'- 


r 


km' 10 kpc 


kT 
(keV) 


Flux2 


C^ (dof)3 


Chandra ObsID 104 
XMM ObsID 0111070201 
Chandra ObsID 3504 
XMM ObsID 0152620101 


1 3+0.4 

1 qO+0.06 

J--'JO_0.04 


12.1 ±0.5 
8.0lgi 
8.4 ±0.2 


0.88 ±0.03 
0.80 ±0.01 
0.81 ± 0.02 

n 7H9+0.008 


gQO+400 

3301^0 

gQO+400 

2901^0 


0.21 ± 0.02 

Q 254+0.008 

0.19 ± 0.01 
245+0-003 


2.2 

1.7 
1.7 
1.5 


1.04 (1844) 
1.15 (1888) 
1.10 (1926) 
1.38 (1857) 



"Normalization of power-law component at 1 keV in units of 10 ^ erg cm ^ s ^ keV 

''Absorbed 0.3-10.0 keV flux in units of 10""^° erg cm~^ s~^ . 

"^The reduced Cash statistic and number of degrees of freedom for the fit. 
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Fig. 2.— The X-ray flux history of 4U 1626-67 from 1977-2003. 
The dashed line is a linear fit; the dot-dash line a logarithmic fit 
to the data. Error bars represent the 1 a confidence intervals. The 
most recent four data points represent the Chandra and XMM- 
Newton observations. 



blackbody component is preferred over a single power- 
law at the 8 a level, and this is what we use for our 
continuum fits. To model the absorption, w e used an up- 
dated version of the tbabs model, tbnew^ (|Wilms et alj 
l2007f ). which includes high- resolution structure for the Ne 
K, Fe L, and O K edges. The best-fit continuum models 
are presented in Table [2l and the combined continuum 
spectra are shown in Figure [1] 

The 0.3-10 keV flux followed a decreasing trend, falling 
from 2.2 x fO^^" erg cm^^ s~^ in 2000 September to 
1.5 X 10-1° erg cm'^ s'^ in 2003 August. Overall, 4U 
1626-67 has decreased in flux since 1977, with no appar- 
ent change following the reversal in accretion torque in 
1990 (see Figured]). 

3.2. Emission Lines 

There are prominent Ne and O emission lines in all 
four spectra (see Figures [1] and [3]) . To explore the line 
characteristics, we added Gaussian features to the con- 

^ http:/ /astro. uni-tuebingen.de/'^wilms/research/tbabs/ 
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Fig. 3. — Emission line regions for all four datasets. The 
hydrogen-like Ne X (left) and O VIII (right) ions are in the top 
panels, and the helium-like triples of Ne IX (left) and O VII 
(right) are in the bottom panels. Each panel contains, from top to 
bottom, Chandra ObsID 104, XMM-Newton ObsID 0111070201, 
Chandra ObsID 3504, and XMM-Newton ObsID 0152620101. 



tinuum models determined from the overall spectral fits. 
All the lines are well- fit by single Gaussian profiles. We 
note that the Ne IX lines are not reliably detected in 
the XMM-Newton data; for these lines, we fixed the line 
width and location to the Chandra values to obtain ap- 
proximations of the line flux. The single- Guassian line 
fits are presented in Table [3] 

The single-Gaussian fits reveal lines which are quite 
broad (« 6000 km s^^ FWHM). For a given hnc, the 
width does not vary significantly over the four observa- 



5 



TABLE 3 
Single-Gaussian Emission Line Fits 



Instrument 


MJD 


AA 

(A) 


V 

(km s"-*-) 


FWHM 
(km s"!) 


Fluxl 


EqW 
(A [eV]) 


C^ (dof)2 


Ne X line at 12.13 A 


Chandra 
XMM-Newton 
Chandra 
XMM-Newton 


51803.6 
52145.1 
52795.1 
52871.2 


0.021 ±0.014 

016+° °^^ 
U.UiD_Q 

O.OlS_o 
019+° °2i 


510 ±340 

'-'='"-880 
^^'J-320 
^"^-540 


6330+""" 

DOOU_ggo 

4800li«oo 
55401^^" 
5540l-«° 


41 s+i-o 

25.011:? 
19.4l|2 

20.411" 


0.29 [24] 
0.21 [17] 
0.20 [17] 
0.21 [18] 


0.97 (115) 
1.22 (16) 
1.04 ( 134) 
2.15 (16) 


Ne IX line at 13.55 A 


Chandra 
XMM-Newton^ 
Chandra 
XMM-Newton^ 


51803.6 
52145.1 
52795.1 
52871.2 


010+'^-"^'' 
U.U1U_Q 032 

u.uiy_o 039 


^^'^-710 
^'^'^-850 


4200+^**"" 

4200 
45QQ+2000 

4500 


13.2l-» 
3 4+4.6 

7 n + 2.6 
'■^-2.5 

'■1-2.1 


0.10 [6.8] 
0.030 [2.0] 
0.079 [5.4] 
0.070 [4.7] 


1.02 (91) 

1.03 (17) 
0.90 ( 102) 

0.81 (17) 


O VIII line at 18.97 A 


Chandra 
XMM-Newton 
Chandra 
XMM-Newton 


51803.6 
52145.1 
52795.1 
52871.2 


U.UDJ_o 079 
U.UZO_o 032 

-0 01fi+° °'5'' 
U.UlD_o 060 

-0 006+° °^^ 
U.UUD_o 021 


1000+^™ 
iUUU_i250 

'^^'^-500 

-260 ± 940 

J^'JU_34o 


(iUU_i7oo 
00±U_ggo 

6000l?«0 
5840l«?° 


601}^ 
SS.llsl 

ZS.»_g 2 

26.ll2;» 


0.56 [19] 
0.49 [17] 
0.40 [14] 
0.40 [14] 


1.76 (45) 
1.43 (60) 
0.99 ( 56) 
1.34 ( 61) 


O VII line at 21.80 A 


Chandra 
XMM-Newton 
Chandra 
XMM-Newton 


51803.6 
52145.1 
52795.1 
52871.2 


-0 050+"-"''* 
U.UOU_o 077 

-0 027+°-"™ 

-0 020+° "^^ 
U.UZU_Q 058 

U.U01_Q 023 


-688+-- 

Qon+960 
oou_i280 
9on+1280 
-^«U_goo 

'^^^-320 


6154lir7 
8900lf~ 

O70Q+3600 
"J" -1400 

52001850 


'7'Y o^h-^b.ii 
' ' ■'^-24.3 
79+18 
'^-13 
3,+16 
■^1-12 


0.87 [23] 
0.96 [25] 
0.50 [13] 
0.86 [23] 


0.98 (22) 
1.08 (50) 
0.91 ( 25) 
1.18 ( 51) 



"The Gaussian normalization in units of 10~^ photons cm~^ s~^. 

''The reduced Cash statistic and number of degrees of freedom for the fit. 

"^Since the XMM data were not able to constrain the Ne IX lines, their positions and widths were 
fixed to the values found in the Chandra data. 



tions. However, the line strengths of Ne X and O VIII 
decrease over the course of the observations (the trend is 
not so clear for Ne IX or O VII). To determine whether 
the line flux is decreasing more rapidly than the overall 
continuum flux, we fixed the line widths to their weighted 
averages and recalculated the line strengths, then fit the 
fiux trends for the Ne X and O VIII lines. Fixing the 
widths did not change the fiux values appreciably, and 
improved the errors on the fiuxes by only a small amount. 
Over the « 3 year span of the observations, from 2000 
September to 2003 August, the Ne X fiux decreased to 
51.6 ± 5.5% and the O VIII fiux to 54.0 ± 8.7% of their 
initial values. This decrease is significantly more than 
the decrease in the 0.3—10.0 keV continuum fiux, which 
only fell to 76.4 ± 6.6% of its initial value. The fact 
that the line strengths decrease faster than the contin- 
uum strength is also refiected in the declining equivalent 
widths of the lines. We find that the equivalent widths of 
Ne X and O VIII decrease to 75.7±7.9% and 74.3±12.3% 
of their initial values, respectively (see Figure IH which 
also includes data from ASCA and BcppoSAX). 

The Chandra Ne X and O VIII lines, as well as the 
XMM-Newton O VII line, show what appear to be 
double-peaked profiles (see Figure [31) ■ As SOI point out, 
the line shapes, combined with the apparent high-density 
environment of the line-formation regions (as indicated 
by the dominance of the intercombination component of 
the He-like triplets; see below), suggest that the lines 
arise in an accretion disk. To derive physically meaning- 
ful values for the line parameters, we fit each line with a 
pair of Gaussians, fixing the width and absolute velocity 



of the red- and blue-shifted components to be equal, thus 
approximating the expected disk-line "two-horn" profile. 
However, there are a number of instances where it is not 
possible to constrain the width and velocity of a partic- 
ular line. For these instances, we fixed the values to the 
weighted averages from the observations where the line 
was well constrainted. This allowed us to evaluate the 
approximate fiuxes of the line components, with the as- 
sumption that the widths and velocities of a given line 
do not change substantially over the course of the obser- 
vations. Where multiple measurements could be made, 
this is true both for the widths from the single-Gaussian 
fits and for the widths and velocities from the double- 
Gaussian fits, so we consider it to be a reasonable as- 
sumption. The Ne IX lines in the XMM-Newton data 
were too poorly constrained to be included in the double- 
Gaussian fits. Table m contains the values obtained from 
these fits. 

Fitting the helium- like triplets of Ne IX and O VII with 
multiple Gaussians allowed us to estimate the strengths 
of the resonance, intercombination, and forbidden lines. 
The Ne IX intercombination line is clearly dominant in 
both Chandra observations (we were unable to fit this 
line in the XMM-Newton observations) . While the inter- 
combination line is the strongest of the O VII triplet as 
well, the resonance line contains significant flux, whereas 
the forbidden lines do not contribute substantially to any 
of the triplets. The lack of forbidden lines suggests that 
we are observing emission from a high-density plasma 
(lie > 10^^ cm~'^), which in turn supports the hypothesis 
that the lines arise in the accretion disk. 
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TABLE 4 
Double-Gaussian Emission Line Fits 













Blueshifted Lines 


Redshifted Lines 








FWHM 


AA 


V 


EqW 


EqW 




Observatory 


MJD 


(km s-l) 


(A) 


(km s-i) 


Flux^ (A [eV]) 


Flux^ (A [eV]) 


(dof)2 



Ne X line at 12.13 A 



Chandra 


51803.6 


36901630 


±0.075±°;°il 


±1840^^™ 


I5.8I3 g 


0.104 


[8.84] 


23.8l3;| 


0.160 [13.36] 


1.02 (246) 


XMM-Newton 


52145.1 




±0.077l°;°l? 


±1910l«° 




0.092 


[7.82] 


14.415-2 


0.120 [9.99] 


1.15 (33) 


Chandra 


52795.1 


3040t^I° 


±0.068 ± 0.007 


±1670+i?E; 


8.2lii 


0.083 


[7.05] 


lo.sll;! 


0.107 [8.94] 


1.10 (281) 


XMM-Newton^ 


52871.2 


3085 


+0 072+°-°l'' 


±1780tg° 


11.3t|« 


0.113 


[9.45] 


9.0l2-j 


0.089 [7.60] 


1.86 (34) 


Ne IX resonance line at 13.45 A 


Chandra'^ 


51803.6 


1300 


±0.055 


±1210 


1+1-' 


0.001 


[0.05] 


1+1-"^ 


0.001 [0.05] 


1.09 (113) 


Chandra 


52795.1 


1300+9^° 


+0 055+°-°l'' 


=ClZlU_4gg 


1 1+1-° 


0.012 


[0.80] 


1 1+1-° 


0.012 [0.80] 


0.86 (120) 


Ne IX intercombination line at 13.55 A 


Chandra^ 


51803.6 


1300 


±0.055 


±1210 


'^•'-'-1.9 


0.037 


[2.52] 


5 0+1-** 
''•'-'-1.9 


0.037 [2.52] 


1.09 (113) 


Chandra 


52795.1 


1300+^^° 


±0.055«;°1^ 


+ 1210+'^*° 

=ClZlU_4gg 




0.021 


[1.46] 


2-311;? 


0.021 [1.46] 


0.86 (120) 


No IX forbidden line at 13.70 A 


Chandra"*^ 


51803.6 


1300 


±0.055 


±1210 


2+1-' 


0.001 


[0.09] 


2+1-' 
'-'•^-0.2 


0.001 [0.09] 


1.09 (113) 


Chandra 


52795.1 


i3oo±;^?° 


±0.055«;°iJ 


ltiZiU_4gg 


0+1-2 


0.0 [0.0] 


0+1-2 


0.0 [0.0] 


0.86 (120) 


O VIII line at 18.97 A 


Chandra^ 


51803.6 


2474 


±0.116 


±1840 


26.31-^ 


0.258 


[8.78] 


21.6li°4« 


0.210 [7.31] 


1.72 (50) 


XMM-Newton 


52145.1 


2450+900 


±0 122+°-°l^ 


±1930 ± 260 


17.6tll 


0.217 


[7.58] 


21.9li^ 


0.275 [9.36] 


1.10 (73) 


Chandra 
XMM-Newton 


52795.1 
52871.2 


2320+1140 

2fion+*™ 

ZDUU_gQQ 


±0 112+°-°2l 
ii:u.iiz_g 023 

+0 n 5+0-012 


±1770l37g 
±1810j:i»° 


13.015-7 


0.179 
0.191 


[6.24] 
[6.66] 


1 •? 7+5.8 
l'^- ' -5.0 

12 4+2° 
1^-^-1.8 


0.191 [6.51] 
0.189 [6.43] 


0.98 (61) 
1.15 (76) 


O VII resonance line at 21.60 A 


Chandra^ 


51803.6 


1750 


±0.104 


±1430 


12.811- 


0.142 


[3.82] 


12.8lllg2 


0.142 [3.82] 


1.16 (22) 


XMM-Newton'^ 


52145.1 


1750 


±0.104 


±1430 


11.7l^i 


0.171 


[4.60] 


11.71^:^ 


0.171 [4.60] 


1.16 (47) 


Chandra^ 


52795.1 


1750 


±0.104 


±1430 




0.048 


[1.30] 


3.ltL? 


0.048 [1.30] 


0.95 (25) 


XMM-Newton 


52871.2 




+0 in4+°-°i'' 


±1430 ± 220 


f\ Q+2.3 
0-^-0.6 


0.122 


[3.27] 




0.122 [3.27] 


0.96 (49) 


O VII intercombination line at 21.80 A 


Chandra^ 


51803.6 


1750 


±0.104 


±1430 


19.8t-- 


0.110 


[2.89] 


I9.8l-;g 


0.110 [2.89] 


1.16 (22) 


XMM-Newton^ 


52145.1 


1750 


±0.104 


±1430 


15. 3+71 


0.100 


[2.63] 


15. 3+71 


0.100 [2.63] 


1.16 (47) 


Chandra^ 


52795.1 


1750 


±0.104 


±1430 


13.6l?-° 


0.159 


[4.19] 


13.6+?-S 


0.159 [4.19] 


0.95 (25) 


XMM-Newton 


52871.2 




+0 104+0016 

ltU.lU4_g gj^g 


±1430 ± 220 


16.611;^ 


0.156 


[4.11] 


16.611;^ 


0.156 [4.11] 


0.96 (49) 


O VII forbidden line at 22.10 A 


Chandra^ 


51803.6 


1750 


±0.104 


±1430 


9 o+iU.U 
^■°-2.8 


0.020 


[0.52] 


r, o + iO.O 

^-°-2.8 


0.020 [0.52] 


1.16 (22) 


XMM-Newton'^ 


52145.1 


1750 


±0.104 


±1430 


5.7l^;| 


0.057 


[1.46] 


5-7l|:l 


0.057 [1.46] 


1.16 (47) 


Chandra^ 


52795.1 


1750 


±0.104 


±1430 




0.0 


0.0] 


o.ol^g;S 


0.0 [0.0] 


0.95 (25) 


XMM-Newton 


52871.2 


1750+^^° 
J-'OU_44o 


+0 104+0-017 

ltU.iU4_g gj^g 


±1430 ± 220 


0+'^-^ 

'J-^-o.o 


0.0 


0.0] 


0+2-^ 


0.0 [0.0] 


0.96 (49) 



"^Thc Gaussian normalization in units of 10^^ photons cm^^ ^-1 

^The reduced Cash statistic and number of degrees of freedom for the fit. 

^The Ne X line -width was not well constrained for this observation, and its value was fixed to the weighted average of the values found in the other observations. 
'^The Ne IX line widths and velocities were not well constrained for this observation and were fixed to the values found in the Chandra observation of MJD 52795.1. 
"^The O "Vlll line width and velocity were not well constrained for this observation and were fixed to the weighted averages of the values found in the other observations 
-^The O -VII line width and velocity were not well constrained for this observation and were fixed to the values found in the XMM-JVewton observation of MJD 52871.2. 



3.3. Photoelectric Absorption and Elemental 
Abundances 

Since 4U 1626-67 sho-ws prominent Ne and O emis- 
sion lines, it is natural to ask -whether this is a reflec- 
tion of an overabundance of heavy elements in the sys- 
tem, or of the particular geometry and conditions in- 
herent to this unique source. One -way to constrain the 
local abundances is by measuring the strength of the ab- 
sorption complexes -within our spectral range. For the 
XMM-Newton and Chandra data this includes the Ne 
K , Fe L, and O K edges. We fit the edges over the 
wavelength ranges 13.8-14.8, 16.5-18.0, and 22.2-25.3 
A, respectively, and used the tbnew model, -which in- 
cludes high-resolution modehng of the edge cross-sections 
(|Wilms et al .1 [20071 : [JueFt et al.|[200l [200l . In order to 



constrain the continuum over the restricted -wavelength 
ranges, we used the same continuum model as -was fit- 
ted to the full spectrum, restricting the model parame- 
ters to vary -within the previously determined 90% con- 
fidence intervals. We present the fit results in Table O 
Although our results suggest a slight overabundance of 
these elements relative to the expected interstellar value 
of Na = (6.2 ±0.7) X 10^° cm~^ (measured directly from 
the Lya line; I Wang fc Chakrabart-yl [200^ . neither the 
Chandra nor the XMM-Newton datasets require a sta- 
tistically significant local contribution to the interstellar 
column depth (see Section [57T] for further discussion). 



4. TIMING ANALYSIS 



7 



Yeor 

1992 1994 1996 1998 2000 2002 2004 



Year 

1980 1985 1990 1995 2000 2005 



60 
50 
40 
30 
20 
10 
2% 
20 
15 
10 
5 
48 
30 
20 
10 
68 

401- 



5 20 
o 



^ <fll 



49000 



50000 51000 52000 
Julian date - 2,400,000.5 



53000 



Fig. 4. — Line equivalent widths from the Uterature and from 
the single-Gaussian fits presented in this paper. Th e ASC A mea- 
surements are plotted with stars l|Angelini et al.l 119951) , Bep- 
poSAX with a triangle l lOwens et al.l 119971) . Chandra with dia- 
monds, and XMM-JVewton with squares; error-bars represent the 
90% confidence intervals. 



TABLE 5 

4U 1626-67 ABSORPTION EDGE FITS 



A 

(A) 


Edge 


Observatory 
(ObsID) 


Total A^H^ 
(10^1 cm-2) 


Implied Nz^ 
{10^"^ cm'2-) 


14.3 


Ne K 


Chandra (104) 
XMM (0111070201) 

Chandra (3504) 
XMM (0152620101) 


i.9t;? 


2.6+0-^ 
l.UO_Q 29 


17.5 


Fe Liii 


Chandra ObsID 104 
XMM (0111070201) 

Chandra (3504) 
XMM (0152620101) 


6+2-^ 
'^•"-0.3 
-1 00+0. 70 

i.3»_Q 



u.a_Q 2 

1 2Q+0.35 


-0.07 
Q 07-1-0.19 
'-'■■"-0.05 

'-'■l'^-0.05 
U.34»_o 


23.3 


K 


Chandra ObsID 104 
XMM (0111070201) 

Chandra (3504) 
XMM (0152620101) 


, 00+0. 42 

1 7-1-0.6 
l-'-1.4 
, 00-1-0. 25 
l-3S_o.o4 


5 9+'''-^ 

^•^-3.4 

6 7+2-^ 
6.8tJ-2 



Fitted value of N^, using tlic tbnew model fWTT 
imediatc range of the absorption structure, restricting the 
main within their 90% confidence intervals (see text). 
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Fig. 5. — Frequency history of 4U 1626—67. Top panel 
shows all data; bottom panel is an expanded view of the Comp- 
ton/BATSE through the Chandra and XMM-JVewton observations, 
along with the BATSE prediction (solid curve with dot-dash er- 
rors). Pulse frequency data from othe r observatories, as well as 
the BATSE prediction, can b e found in IChakrabartv et al.l (119971 ') 
and lChakrabartv et all 1 120011) . 



All of the Chandra and XMM-Newton observations of 
4U 1626-67 have high enough time resolution to deter- 
mine the pulse period, and the XMM-Newton observa- 
tions also allowed us to create energy-resolved pulse pro- 
files. Prior to performing the timing analysis, we cor- 
rected the photon arrival times to the location of the 
solar system barycenter. 

To determine the pulse period, we used an oversampled 
fast Fourier transform, with the frequency and frequency 
error determined by the following formulae ([Chakrabartvl 
fT996tlMiddleditchlil976l:rMiddleditch fc Nelsonlll976lh 



47r2er 



P, 



Pk 



(1) 



0.014 r 

0.012 j- 

7 0.010 r 



10.0- M.0 kev 
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in 

I 0.06 
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0.20 : 
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Fig. 6.— Folded pulse profiles of 4U 1626-67 from the XMM- 
Newton PN, ObsID 0111070201. 
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Fig. 7.— Folded pulse profiles of 4U 1626-67 from the XMM- 
Newton PN, ObsID 0152620101. 



2ttT 



(2) 



where P is the normahzed power, the fc*'' bin con- 
tains the peak power, T is the length of the observa- 



tion, and 



e is the oversampUng factor (for our 



analysis we used n — 6). The results are presented in 
Table IH and Figure [5] shows a plot of the frequencies 
as well as the prediction curve and associated errors ex- 
trapolat ed from previous Compt on/BATSE monitoring 
results ([Chakrabartv et al]|1997f) . As can be seen from 
the plot, the pulsar is spinning down faster than was pre- 
dicted, which is significant with respect to the stability 
of the timing during the BATSE era, but not unexpected 
considering the variable spin history of the source. We 
note that the QPO at 0.048 Hz is present in eac h of the 
observations, as has been seen previously fShino da et al.l 
[l990t lChakrabartvlfT99l . 

The XMM PN data are ideal for creating energy- 
resolved pulse profiles. We made four energy cuts and 
divided the pulse period into 40 bins to create the profiles 
presented in Figures [6] and [71 The profiles have some sig- 
nifi cant differences from what was seen previously (see, 
e.g.lPravdo et al.lll979l:lKii et al.lll98llMcCHntock et all 
ll980HLevine et al.lll988[ ). In the earlier observations, the 



pulse profile shows a dip at energies < 2 keV, then the 
emergence of two prominent peaks which flank this dip 
as the energy increases, reaching a maximum at ~ 5 keV. 
By « 13 keV, the dip has returned, and becomes both 
broader and deeper as the energy increases. We find 
that the current pulse profiles lack the double-peak fea- 
ture that was seen before, and show only a dip which 
broadens at lower and higher energies. The profile is not 
entirely symmetric; there also appears to be a secondary 
dip at phase « 0.85. Insofar as the pulse profile reflects 
the geometric and radiative properties of the accretion 
colum n and neutron star hot spot (see, e.g. iKu et all 
Il986f ). changes in the pulse profile suggest that there 
have been fundamental changes in these regions of the 
system. While we do not speculate here as to the exact 
nature of these changes, we note that the accretion ge- 
ometry may have changed at the time of torque reversal, 
affecting the accretion fiow at the neutron star surface. 
Therefore, the changes that are seen in the pulse profile 
may be related to the change in the continuum spectrum 
of 4U 1626—67 t hat was observed aro und the time of the 
torque reversal (|Angelini et al.l ll9951. 

We also note that there were no fiaring events seen in 
any of our observations, in contrast to what has been pre- 
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TABLE 6 
Pulse period of 4U 1626-67 







Pulse Period 




Date 


MJD 


(s) 


Observatory 


2000 Sep 16 


51803.6 


7.6726(2) 


Chandra 


2001 Aug 24 


52145.1 


7.6736(2) 


XMM-Newton 


2003 Jun 5 


52795.1 


7.67514(5) 


Chandra 


2003 Aug 20 


52871.2 


7.67544(6) 


XMM-Newton 



viously observed, when 4U 1626—67 was seen to flare dra- 
matically in b oth X-ray and o p tical data on timescales 
of pa 1(300 s (iJoss et all llQTSl: iMcClintock eTaH [l980l: 
iLi et al.lllQSOf ) . The cessation of flaring activity may have 
occurred at the same time as the torque reversal, though 
it is not clear what physical mechanism was responsible 
for the flaring or why it would be correlated with the 
change in sign of P (|Chakrabartv et al.l[200l . 

5. DISCUSSION 
5.1. Absorption Edge Measurements 

Four known or candidate ultracompact binary systems 
have X-ray spectra which show evidence of high Ne/0 ra- 
tios (4U 614-^091, 2S 0918-549, 4U 1543-624, and 4U 
1850-087: Uuett et al.ll200ll : lJuett fc Chakrabartvl[2003l) . 
though none of these systems have the X-ray emission 
line complexes seen in 4U 1626-67. These unusually 
strong lines, which arise from highly ionized species of 
Ne and O, suggest that the donor must contain a sig- 
nificant quantity of these elements, and we have previ- 
ously suggested that the donor may be the chemically 
fractionated core of a C-O-Ne or 0-Ne-Mg white dwarf 
(SOI). However, since the strength of the emission lines 
is very dependent on local plasma parameters, it is dif- 
ficult to use them to determine elemental abundances. 
Alternatively, a way of determining local abundances is 
to measure the depths of the neutral absorption edges 
and to compare these with the value of the interstellar 
absorption along the line of sight (for example, as deter- 
mined by the broader continuum fit). Any absorption 
in excess of what is expected from standard ISM abun- 
dances could be taken as an indication of material local 
to the system. 

The Ne K, Fe L, and O K absorption edges do not 
appear strong in any of our observations, including the 
first Chandra observation (ObsID 104) that we analyzed 
in SOI. Therefore, we are no longer able to argue for 
a significant overabundance of O or Ne in cool circum- 
stellar material as we suggested in SOI, at least not on 
the grounds of the absorption edge analysis (see also Sec- 
tion [3?3l)- Perhaps these elements are, in fact, overabun- 
dant in the disk (and thus the donor), but not much ma- 
terial is expelled from the system, or subsequently cools 
to form an absorbing medium. Our results agree with 
a recent study o f the optical spect ra of 4U 0614-1-091 
and 4U 1626-67 (| Werner et al.ll2006l ). These spectra did 
not reveal any Ne emission features, only lines from C 
and O, which allowed them to place a limit on the Ne/0 
ratio of « 0.2. Furthermore, 4U 0614+091 showed evi- 
dence of a higher C/0 ratio than 4U 1626-67, suggest- 
ing that its donor is actually more evolved. If this is the 



case, than the most significant difference between the two 
systems — the high magnetic field of 4U 1626-67 — could 
be the reason that 4U 1626-67 has prominent emission 
line features whereas the other ultracompact sources do 
not. 

5.2. Emission Line Characteristics and Physical 
Implications 

The emission lines have double-peaked profiles, and the 
line ratios imply that they are formed in a high-density 
environment, which suggests that they arise somewhere 
in the accretion disk (see also SOI). If they are, in fact, 
disk features, then the measured velocities correspond to 
the actual disk velocities by the relation 



i^obs = f sin i = Y ^sinz. (3) 

The maximum disk velocity will occur at the in- 
ner edge of the accretion disk, which is truncated at 
the co-rotation radius of the neutron star at « 6.5 x 
10^ cm. The Keplerian velocity at this radius is ~ 
5400 km s~^, whereas the line velocities are measured 
to be ~ 1700 km s~^ (note that this is the velocity at 
the center of the Gaussian in the double-Gaussian fits, 
which is therefore an underestimate of the absolute max- 
imum velocity). If we take this to be the velocity at the 
inner disk edge, we are able to place a constraint on the 
inclination angle of the system: sini > 1700/5400 ~ 0.31 
(i > 22°). We may combine this with the previously de- 
rived upper limit on the projected semi-major axis, given 
an orbital period of 42 mi n and limits on the timing noise 
(jChakrabartv et al.lll997l ). to find 3 < Ox sini < 8 It-ms. 

The line ratio measurements are indicators of the 
plasma conditions in the line formation regions. Al- 
though we are not able to derive robust limits for the 
standard diagnostics, we note that the lack of detectable 
forbidden lines implies that the region is relatively high- 
density (rio ?J 10^^ cm"'^). Furthermore, the relatively 
low value of the intercombination line (z) with respect 
to the resonance line (r), r/i « 2, suggests that the 
temperature of the line -formation region is > 10^ K 
(|Porquet fc Dubaull2000f l. As was discussed in SOI, this 
temperature is characteristic of the highly ionized, opti- 
cally thin outer layers of the accretion disk. 

5.3. Long-term Flux Evolution 

X-ray observations of 4U 1626-67 show that it has been 
decreasing in flux since 1977. Since this decrease appears 
to be bolometric — there is no sign of variable absorption, 
and the continuum spectrum is known to have undergone 
only one major change — we may use the X-ray flux to 
help determine other characteristics of the system. The 
long-term trend can be described by either an exponen- 
tial decay or a linear decrease (see Figure If we take 
it to be logarithmic, and integrate over the duration of 
the outburst, the fluence is 0.927 erg cm^ (giving a total 
energy of 1.1 x 10^** (d/1 kpc)^ erg). Converting this to 
an average accretion rate, we find that 

Afavc ~ 2.2 X 10-1" At;/ dl^^ Mq yr-i, (4) 

where Atyi- is the time between outbursts in years and 
c?kpc is the distance to the source in kpc. If we assume 
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that the source is persistent and take « 30 yr, we 
find that < 2 kpc, which is less than the minimum 
distance of 3 kpc derived from considera tions of the ac - 
cretion torque during spin-up (fChakraba rtv et aLlll997[ ) . 
If we assume a distance of > 3 kpc, and take the long- 
term average accretion rate to be 3 x 10^"'^^ Mq yr~^ (the 
value from gravitational-radiation- driven orbital evolu - 
tion with a white dwarf donor; see lChakrabartvlfT99l . 
we find that At > 70 yr, significantly longer than the 
known duration of the current outburst. 

We ma y also combine current theory on dis k stability 
(see, e.g. iKing et"all 119961 : Ivan Paradii3ll996f ) with the 
flux data for 4U 1626-67 to obtain constraints on its 
distance and possible outburst duration. If the effective 
disk temperature is that of an X-ray heated accretion 
disk, then 



4 _ riMc^{l-l3) fdH 



H 



dr 



(5) 



Following [King et all (|1996[ ). we take 77 0.11, /3 = 
0.9, assume that H/r ~ 0.2 is constant, and that H oc 
r", where n = 9/8 for shallow heating, which we would 
expect to dominate in this system. The disk stability 



criterion is that the te mperature at the outermost disk 
radius remain > 6500 K ()Smaklll983h . Setting r — rant — 
2 X 10"'^° cm, we derive a critical accretion rate of M = 
3.6 X 10"^^ Mq yr-^ Note that this is lower than Mq^, 
Assuming that the accretion rate is directly proportional 
to the flux, and that it remained above the critical rate 
throughout the current outburst, we find the distance to 
4U 1626-67 to be > 0.6 kpc. Furthermore, if we assume a 
distance range of 3-13 kpc and an exponential (or linear) 
decline in flux commensurate with the current trend, we 
flnd that M < Merit in the range 2018-2020 (or 2008 
June-Sept). We therefore expect 4U 1626-67 to become 
quiescent within 2-15 years. 
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